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1 Introduction 



Solving matrix models using the loop equation has more than decade long history. An 
algorithm for finding genus expansion solutions to the Hermitian one-matrix model was 
elaborated (see pQ and references therein) in the case where the limiting eigenvalue distri- 
bution spans a single interval on the real line in the large N limit (in what follows, N is 
the dimension of the Hermitian matrix). Such a situation is commonly called the one-cut 
solution. Results for multi-cut solutions were few. The problem was that it is not too 
difficult to find a solution for the loop equation, i.e., the genus expansion of the so-called 
loop mean — this problem was solved by Akemann [2|. The universal critical behavior of 
the corresponding correlators was considered in However, a difficult problem is to 
integrate the obtained solution in order to produce higher genus contributions to the free 
energy. Akemann has succeeded in finding such a genus-one contribution in the two-cut 
case only. 

Recently, it was realized that multi-cut solutions to matrix models is the powerful 
tool for studying broken vacua structures of M = 2 super Yang-Mills theory of Seiberg 
and Witten Simultaneously, a new relation between the superpotentials of M = 1 
super symmetric gauge theories in four dimensions and free energies of matrix models in 
the planar limit was proposed [HllSl- The superpotentials in some M = 1 four- dimensional 
Yang-Mills theories has turned out to be expressed through a single holomorphic function 
j^], which was later identified with the planar limit of the free energy of the multi-support 
solutions to matrix models j^j. 

On the other hand, there are strong indications that such an analogy goes beyond the 
tree level and just the large- N limit of the matrix model partition function. The relation 
between gauge theories and matrix model solutions were checked on the one-loop level 
for the solution with two cuts and a cubic matrix model potential [7j. A comparison in 
general case must involve the structure of the general matrix model multi-cut solution in 

1 E-mail: chckhov@mi.ras.ru. 



1 



genus one approximation, which has been missed as yet. Kostov [S] had provided argu- 
ments in favor of using the conformal field theory (CFT) methods for calculating such the 
corrections in the general case following the way proposed by Moore pj; this method has 
been improved by Dijkgraaf, Sinkovics, and Temurhan ^U] who used algebro-geometrical 
and CFT arguments for proposing the formula for genus one correction and comparing 
it with the gauge theory answer. Note, however, that the very applicability of Kostov's 
star operator (dressing) method needs to be justified by direct calculations using the loop 
equation alone and, second, the algebro-geometrical methods need to be improved in a 
way to allow calculating contributions from "non-geometrical" part of the action — from 
the so-called "moments" of the matrix model potential (such contributions always arise if 
the power of the polynomial that is the derivative of the potential exceeds the number of 
cuts in play; then, double points arise and the genus of the arising hyperelliptic surface 
can be much lesser than the highest degree of the potential. In this situation, however, the 
superpotential satisfying the Whitham hierarchy equations and the WDVV equations is 
still attained by the planar (large- N) limit of the free energy of the corresponding matrix 
model [HI, Cl- 
in this paper, we derive the subleading (genus one) expression for the free energy of 
the generalized matrix model multi-cut solution directly from the loop equation, without 
referring to geometry or to the CFT reasoning thus supporting the calculations in jH] 
and [TP] . In Sees. 2, 3, we closely follows the method of [I], previewing the iterative 
procedure for finding the loop mean in the genus expansion; the answer for the genus-one 
contribution W\(jp) to the loop mean is presented in Sec. 4. We integrate the obtained 
answer in Sec. 4 thus obtaining the complete expression for the genus-one part F\ of 
the free energy. As a byproduct, we obtain interesting combinatorial identities on the 
quantities entering Wi(p). 

2 Definitions. Loop equation 

The partition function of the Hermitian one-matrix model is 

Z[N,{ti}] = e N2p = e^7=o N2 - 29F ^ = J dSe-^W (21) 

where the integration is over Hermitian N x N matrices $. The matrix potential, which 
we assume to be polynomial of degree (m + 1), is given by the power series 

m+1 j 

V(x) = ~ xJ - (2-2) 
j= i J 

Deriving the loop equation implies introducing all coupling constants tj into play; a specific 
potential of finite order can be inspected by setting the extra couplings to zero in the final 
result. The coupling constants tj are then the sources for polynomial expectation values 
with the means defined customarily as 

(Q{<l>)) = ^fd*Q(*)e- N * v <». 
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Introducing the loop insertion operator 

>> = -£-^, P-3) 
the one- loop mean W(p) can thus be obtained from the free energy F 

N ^ p k+1 N \ p-<j>/ dV yF! p K J 

Analogously, all the multi-loop correlators can be derived by applying J^O) to F (or to 
W{p)) 

W( Pl ,...,p n ) = iV«- 2 (tr^----tr- ' ^ 



pi - <p p 



n 



conn 



Here conn refers to the connected part. Because the loop correlators and the free energy 
have the same genus expansion, 

OO 1 

W( Pl , ...,p n ) = Y^ j^W g ( Pl , ...,p n ), (2.6) 

Eq. ()2.5|) is valid for each genus g > separately. We assume W{p) to be analytic at 
infinity with the asymptotic behavior 

lim W(p) ~ - , (2.7) 

p— too P 

i.e., 

W g {p)\ p ^ oo = U gfi + 0(l/p 2 ). (2.8) 
This means that higher W g (p) must be total derivatives: 

W g (p) = ^(p)F g , g>l. (2.9) 

The loop equation is the exact equation satisfied by the loop mean Deriving the 
loop equation is irrelevant to the multiple cut structure and relies only on the invariance 
of the partition function under a field redefinition $^$ + e/(p — $). The only effect is 
that the contour C-d encircles now a number n of disjoint intervals on the real line: 

{ ^¥^W(w) = {W{p)f + ^-^-(p)W{p\ p?V, (2.10) 
Jc-d 2m p — u JS z dv 

where V'(u) = ^AjUJ^" 1 . We assume that the density Ptv(A) = ^(ZI^^A — ^)) °f 
eigenvalues of the matrices $ has a compact support T> as N — > oo. We let T> to comprise 
n disjoint intervals: 



V = \J\pl>2i-l, fJ>2i], A*l < fJ-2 < ■ ■ ■ < \i-2n- (2- 11 ) 
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Presenting W(p) in terms of pjv(A), W(p) = J e?A PJV _y , we find that as iV — > oo, W(p) is 
analytic outside n cuts on the real axis. 

Inserting genus expansion ()2.6|) into loop equation ()2.1Ujh we obtain 



c B 27ii p — u 



W {u) = (W (p)) 2 (2.12) 



for genus zero and 



{t - 2W (p))W g (p) = J2 W gl {p)W g - gl {p) + ^(p)WV-i(p), 5 > 1, (2-13) 
for higher genera, where /C is a linear integral operator 

*/<p) - / (2-14) 
JCt, 2m p — uj 

Given Wo(p), one can then determine W g (p) for g > 1 iteratively genus by genus provided 
the operator (/C — 2K / o(p)) can be inverted uniquely. 

We can solve Eq. (j2.12|) for the planar solution Wq(p) as follows. Deforming the 



contour in Eq. (j2.12j) to infinity, we obtain 

(W (p)) 2 = V'(p)W (p) + I ^^-W (uj). (2.15) 

J Coo 2m p — uj 

The solution is then formally 

Wo(p) = \v'{p)-\-J(V'( r )f+iP{p), 

P(p) = f pLYMwM, (2.16) 
J Coo 2m p — uj 

where the minus sign is chosen in order to fulfill the asymptotic Eq. (|2.8|) and P(p) is a 
polynomial of power m — 1. If Wo(p) has n cuts in the complex plane, we propose the 
ansatz for the square root, 

W (p) = ^(V'(p)-y(p)), (2.17) 

where 



y(p) = M(p)y(p), and y(p) = \/Y[. =1 (p - fa) (2.18) 

(the notation comes from JI]), and M(p) is assumed to be a polynomial of degree m — n, 
which still has to be determined. By convention, we set y{p)\ p ^oo ~ P n , and M{p) is then 

M( P ) = l (2 . 19) 
J Coo 2m [uj - p)y{uj) 
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Inserting this solution in Eq. (J2.17j) and deforming the contour back, we obtain the planar 
one-loop correlator with an n-cut structure, 

2 jc-d 2m p — uj y(u) 
Then, the planar eigenvalue density p(A) = linijv^oo Pn (A) has the form 

p(A) = lim (w (A - ie) - W (X + ie)) = lmy{\), XeV. (2.21) 

2m e->0 \ 7 Z7T 

We still must determine the positions of branching points p a , a = 1, . . . , 2n. Here, we 
have two sorts of restrictions. The first set comes from asymptotic conditions ()2.8)1 to be 
satisfied by ansatz ()2.2()j) . from which we have 

1 / duu k V'(u) , n , nnn . 
$k,s = -f — - ( \ \ k = 0,...,n. 2.22 

2 JCt, 2m y[oJ) 

These conditions provide n+1 equations for 2n constants fi a (which suffices only if n = 1). 

Another bunch of relations cannot be read from the analytical structure alone. Intro- 
ducing the occupation numbers 

s '=lM M< ^ )s lM y< ^ (2 - 23) 

where Ai, i = 1, ... ,n — 1 is the basis of A-cycles on the hyperelliptic Riemann surface 
y 2 = riali^ — Pa) (we may conveniently choose them to be the first n — 1 cuts), we can 
consider these numbers as independent parameters of the theory in the Dijkgraaf-Vafa 
setting, thus imposing the restrictions 

-£-(p)Si = 0, z = l,...,n-l. (2.24) 
aV 

(Adopting a quite opposite standpoint, which is intrinsic for an original approach to 
the matrix models, see |2j, we must look for the genuine minimum of the matrix 

model action with respect to variables Si as well as with respect to the times tj. This 
means that partial derivatives of Fq with respect to Si must vanish. One can obtain 
(see [Hj, [Hj), that these derivatives are differences of chemical potential on disjoint cuts 
and those are equal to integrals over dual 5-cycles on the corresponding hyperelliptic 
Riemann surface, that is, this set of conditions reads 

-^{p)j B duy{uj) = 0. (2.25) 

Note that, as we shall see below, from the technical viewpoint, it does not matter which 
cycle set, A- or 5-cycles, we take as the basic set. That is, the higher genus corrections 
become scheme dependent: choosing determining conditions ()2.24|) or (|2.25J) . we obtain 
different expressions for the free energy in the the genus one approximation. These ex- 
pressions are however of the same type; the difference will be exactly in the choice of the 
corresponding basis cycles.) 
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3 The iterative procedure 



3.1 Moments of the model 

We can determine higher genus contributions iteratively by inverting genus expanded loop 
equation (|2.13jl . From W g (p), all multi-loop correlators of the same genus can be obtained 
then simply by applying the loop insertion operator ^(p) to it (see Eqs. (|2.5jl and (|2.4|) ). 
As in the one-cut solution pQ, changing variables from coupling constants to moments 
allows expressing higher genus correlators nonperturbatively in the coupling constants tj, 
and these correlators turn out to depend only on a finite number of the moments 

M « ] = L^'^^^ 1 fc = 1 ' 2 '---' « = 1 "..A (3.1) 

with 

0L fc) M = 1 , fc (o) M, (O) H = ^. (3.2) 

Note that expression f)2. 19j) implies 

M« = M(/0- (3-3) 

3.2 Determination of the basis 

We now introduce a basis for the operator acting on W g (p) in Eq. (j2.13)l . 

{K-2Wo( P ))x i a ) (p) = f fc = l,2,..., a = l,...,2n. (3.4) 

Given the r.h.s. of Eq. 1)2.13)1 to be a fractional rational function of p having poles at the 
H a only, W g (p) must have the following structure: 

3g-l 2n 

W g (p)=J2T,^ k ix i a k) (p), g>i. (3.5) 

k=l a=l 

Here A™ are complicated functions of fip and the moments As the order of the 

highest pole in W g (p) is insensitive to a multi-cut structure, W g (p) will depend on at most 
2n(3g — 1) moments, just like the one-cut solution case pp. 
A set of basis functions fulfilling Eq. (j3.4j) is defined by 



^)( p ) = ^(^)(p)- k ^Mt r +^x { : ) (p)], k=l,2,..., a = l,...,2n, (3.6) 

r=l / 



which can be proved by induction. However, this definition is not unique, as the kernel 
of (/C — 2Wo(p)) is not empty, 

Ker(/C - 2W (p)) = Span{pV (0) (p); i = 0,...,n}. (3.7) 
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This follows from boundary conditions (J2.22j) . Now, because of the asymptotic behavior 
()2.8|) . only terms ~ C(pO & t large p may be added to W g (p) for g > 1. This requirement 
reduces the set of zero modes to be 

allowed zero modes : {p l <fi (p)', / = 0, . . . , n — 2} . (3.8) 

In the general case, any linear combination of n — 1 functions can be added to W g (p). 
But, besides its asymptotic behavior, W g (p) must fulfill Eq. 1)2.9)1 . which can be used 
to fix the basis uniquely. The p-dependence of W g (p) must be completely absorbed into 
derivatives with respect to jb(p). This is possible only if the basis functions Xa ip) can ^ e 

(fc) 

expressed completely in terms of ^f(p) and dh ^y fo) as functions of p. So, Xo^(p) mus t be 
redefined in Eq. 1)3.6)1 in accordance with this demand, which fixes completely their zero 
mode components. The derivatives of branching points p, a and the moments M^f 1 can 
be obtained by applying the loop insertion operator ^f(p) to definitions (j3.1j) and using 
restrictions (j2.24j) taking into account that 

Using the identity 

= -<TW' (3 - 10) 

the result for the moments reads 

^fip) = (A: + l/2)(il4 fc+1) ^(p)-^ +1) (p)) 

+ - y y - u®( P ) - 

2n 



2 fri (up - n a ) k \ p dV 



a = l,...,2n , k = 1,2,.... (3.11) 
The quantities ^f(p) are given by the solution of the following set of linear equations, 



= E^^W-^CP))^ J = l,...,n-1, 

if = / rfAf^- , (3.12) 

Jj4j A — /X a 

when applying ^tf(p) to ()2.22j) and using ()2.24|1 . While the first relations can be attained 
easily, deriving the latter can be performed in a tricky way first proposed by Akemann. 
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Indeed, 

d 1 2n 1 du d 

= - 2 E "W^z^w W + »(^WMM (3-13) 

Calculating the second term, we find from (j2.19j) that 



d , , r d£ 



—MM{u) = i 
dv JCoo.p 2.m 



V'{£) fl * 1 dfi 



and the first term can be evaluated by the residue at infinity replacing V'(£,) by M(£)y(£) 
in this limit, while the second term is given by the residue at £ = p (as the integration 
contour must encircle the infinity and p). The first term then becomes 

ifi MM-My dn a/ 

;L ^7W> 3 - 14 

and it partially cancels the first term in 1)3.13)1 . the difference is just the term proportional 
to M(/i a ) = M£\ while evaluating the second term upon substituting it into the integral 
in 1)2.24)1 needs partial integration: 

» i i.-L.M = if i.-^ti— -YM 

dp JAj p-uy{p) 2 J a, p-u ^\uj - fi a p-ii a Jy{p) 

1 2n 

= ^EM?^). ( 3 - 15 ) 

Z a=l 

which together with (jHUBj) and ()3~T4]l gives ()3~T2l . 

Linear system of equations ()3.12jl implies that the solution always takes the form 

du n ~ 2 
M«-^(p) = 0«(p) + £ C a , lP l ^ \p), a = 1, . . . , 2n. (3.16) 

The quantities £ a ,z, a = 1, . . . , 2n, Z = 0, . . . , n — 2, are then determined by the following 
set of equations: 

2n n-2 2n fe-1 

= EE^^f'-EE^'V+s^- 1 , fc = o,...,n, 

a=l 2=0 a=l Z=0 

2n n-2 

= EE £ «,i V. j = l,...,n-l. (3.17) 

We can redefine basis (13.6)) in terms of total derivatives. Solving Eq. 1)3.11)1 for 0^ fc+1 ^(p) 
and expressing there the p-dependence through ^7( P )-terms plus zero modes, we can obtain 
a unique basis inductively from Eq. ()3.6)) properly subtracting the corresponding zero 
modes. We then have 



X { a k) ( P ) = -4y {€\P) ^_ -EMi k - r+1 h^(p)), a = l,...,2n, n = l,2, 



Mi 1 ) 



w -part 



(3.18) 
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Few first basis functions are 



Xa(p) = jy(p)' a = l,...,2n, 



9 (7 1 r\ 

xt\p) = -^ W (p)log\M^\--j:—(p)log\ f i a -^\. (3.19) 



4 Calculations in genus one 
4.1 Genus one partition function 

Having determined the basis, the loop equation can now be inverted step by step in genus. 
For genus g — 1, Eq. ()2.13|) reads 

(/C - 2W (p))W 1 (p) = A(p )Wb (p). (4.1) 

Given Wo(p) (|2.20|) and the loop insertion operator f)3.9j) . the r.h.s. becomes 



a<f3 

+ - 1 V 1 M'^fp) 
4 0(o)(p)^ lP -/i Q ^ 

1^1 1 ^ 1 /1 1 \ 

lQ ^l(p-^a) 2 8 ^ H a - fJL/3 \P - ft a V ~ / 
^ a=l 1=0 V h La 



Here we took into account that regular parts coming from _f , i = 1, . . . , n — 2, vanish 

due to Eq. (|3.17|) in order for Wo(p,p) = -^{p)W (p) to satisfy the correct asymptotic 
behavior, and we can just replace p l by fi l a in numerators of such expressions. The result 
for the one-loop correlator of genus one with n cuts can now be easily obtained using Eq. 

(EDI), 

i 2n i i 

Wk(P) = ^E^GO"* E (X^(P)-X?(P) 

ib «=1 8 l<a</3<2n A*a ~ #9 

i 2n n— 2 

+ I EE^xi 1) (p) 

4 a=l Z=0 
-i 2n i i 

i 2n n— 2 

+iEE^XxL 1) W 

4 a=l Z=0 




3dV KFJ 61 a 1 3 fri dV yFI 



ct,f3 
a<f3 

2ll 



=1 l^a — 1^/3 \ dV dV 



+lE^W^), (4-3) 

where the polynomials C a {p) are, by definition, Ya=o £-a,W l ■ 

We now integrate ()4.3|) in order to obtain F\. While integrating first two terms is easy, 
the term with zero modes is difficult. Even before presenting the result, we need more 
notation. Let us introduce the quantities 

^ i = 1 ' 2 '- (4 ' 4) 

(part of them constitute the ^.-matrix jHj, ^Oj) and the (polynomial) basis of holomorphic 
differentials on the hyperelliptic Riemann surface 

n-1 

H k {\) = Y^H^X 1 ' 1 , k = l,...,n-l, (4.5) 



i=i 



such that 



Ai y{\) 



-d\ = 5k j , (4.6) 



whence 



n— 1 



YsQi,i H i,h = 8j,k for j, = l,...,n- 1. (4.7) 

The following lemma, to be proven in the next subsection, establishes the relation 
between zero modes C a> i and Qj ti . 

Lemma 1 For quantities C a; i determined by system \3.1 r l\j and the quantities \4-4\l an d 
\4-ty > we have the relation 

n—2 n—1 r rr / \ \ 

E^X = -E/ - m ,f J , dX, a = l,...,2n. (4.8) 
l^o Pi J A 3 y(X){X- ii a ) 

It follows from Lemma 1 that the last term in ()4.3j) is just 

2 d/J, a V i,j=l,...,n-l / CfK 

i.e., 

*i = log f ft M« • A 4 ■ (. . det Q^) 12 ) , (4.9) 

where A = ]li<o</3<2n(/ i a — Hp) is the Vandermonde determinant. This is our final answer 
for the genus one partition function. 
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4.2 Proof of the lemma 

Note first that all the quantities Qj t i, although being dependent only on 2n moduli /i a are 
nevertheless algebraically independent for all j and i < 2n (that is, no rational functions 
of n a and Qj^ with i < 2n are identically zero at all /x a ). From the combinatorial point 
of view, we can therefore consider all such Q 3 - ti as independent quantities. We also prove 
relation ()4.8|) not for C a j themselves but for special sums 

2n 

R kti =Y,l4£*,i- (4-10) 

a=l 

If we calculate all Rk,i with k — 0, . . . , 2n — 1, we automatically find all L a ^ through the 
inverse Vandermonde transformation. We can prove ()4.8|) applying the summation with 
/x^, to both its parts. Obtaining identities at all k is equivalent to proving the lemma. 
For shortening notations, it is convenient to introduce the quantities 

n-l 

T ltm =Y,QTjQi,m, 1 = 1,..., n-l, m = l,2,.... (4.11) 

i=i 

Obviously, 7j jm = 5^ m for 1 < m < n — 1. 
In the l.h.s. of ()4.8|) . we therefore have 

2n n-2 

E^«W = E^. ( 4 - 12 ) 

a=l i=0 

while in the r.h.s., we find the expression 

2n n—1 . f\—li—\ k 

Z^xj^jAi y{x)(x-fi a ) 

2n n—1 . Y^ n_1 ( i ft 2n n— 1 /• y n_ 1 1 

_ _ \ -> \ -> / Z^j=l iPo ' _ \ ^ \ " , V £^i=l ^i,] 

a=lj=1 JA 3 y(x)(x-fi a ) a= i j=l JA 3 y(x)(x-fi a ) 

fc-ln-1 / 2n \ n-l 

= EE E /4~ s_1 - 2 E(» + fc - i)W-i> (4.i3) 

s=0 i=l \o=l / i=l 

where we have integrated by part in the second term. Thus, we must prove the coincidence 
of gZEH) and (jHH for all jfe. 

First, it is straightforward to find from the first relations in (|3.17|) the expressions for 
Rk,i with k = 0, . . . , n: 

i? M - 0{k - 1 - z) E - 2^4-^-1,0, fc = 0, 1, . . . , n, (4.14) 

a 

where 8(n) = 1 for n > and zero otherwise. The complementary set of relations follows 
from the second conditions in (|3.17|) . Let us introduce the standard symmetric functions 

Sy m fc = E VhVh ■■■^ k<2n, 

l<il<«2<-..<«fe<2n 

Sym(fi a ) k = E ViiViz • - fJ>i k , k<2n-l 
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with the obvious recurrence relations 

Sym(/4 Q ) fc = Sym fc - fi a Sym(^ a ) fc _ 1 , 

Sym = Sym(^ Q ) = 1, Sym_ 1 = Sym(^ a )_i = 0. (4.15) 

In particular, we have 

Sym(/^ Q ) fc = Sym fc - /i a Sym^ + . . . + (-fi a ) s Sym fc _ s + (-/i Q ) s+1 Sym(^ a ) fc _ s _i 
for s = 0,...,k, (4.16) 

and 

2n 

J2 Sym(^) s = (2n - s) Sym s , < s < 2n. (4.17) 

a=l 

We use the transformation 

il(r) 1 2n 1 2n_1 

= ^n^-M = 4rE^ 1 (-i) i sym(^ 

y{x)^ a y{x) f^ 

j 2n-l / 



V E^ n " ? "VL- s (-l) s Sym s (4.18) 

Z=0 s=0 

to obtain the set of constraints 

2n-2 2n-fc-l 

E E Q,,2„- fc -^Sym fc (-l) fc = 0. (4.19) 

fc=0 s=l 

Using the above formulas and contracting the result with Qpj, p — 1, ■ • • , n — 1, we obtain 
after some algebra using (|4.14|) . (|4.16jl . and (|4.17j) 

2n 



= E^r(-ir n Sym(^) n ^ 1 + (p- ? -l)(-ir 4 - 1 Sym 2 _ 1 _ p _. i 

n— p— 1 

+ £ i4 +s+lii (-l)^ +s + 1 Sym n _ 1 _ p _ s 

s=0 
2n-l 

+ £ T p , r .(-l)^ +1 (r - 1 - i) Sym^^, p = 1, . . . ,n - 1. (4.20) 



We can now find all the remaining quantities R n+ i +s ^ for s = 0, . . . , n — 2 by inverting the 
lower triangular matrix standing by these quantities in the second row of (|4.20|) . Namely, 
we introduce the quantities q r , r = 0, 1, . . ., using the determining relation 



p 

r=0 



J2q r Sym p _ r (-l) p - r = 5 pfi for p = 0,l,.... (4.21) 



The first few q r are q = 1, q%— Sym 1; q 2 = (Syuq) 2 — Sym 2 , etc. Except (I4.21|) . we 
need another combinatorial identity on q r , which follows from (J4.21|) and (|4.15j) . Let us 
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multiply (gZEED by Sy mi = £«=i Mq : 

P 2n p 

5 p,o = E(s + 1) Sym s+1 (-) s g p _ s + E /4 £ Sym(/4 Q ) s _i(-) s g p _ s 

s=0 a=l s=0 

p 2n p+1 tp+X—q 

= E( s + 1) Sym s+1 (-) s g p „ s - EE^l E s y m s(-) s <?P+i-^ 

s=0 a=l <j=2 \ s=0 

p 2n 

s=0 a=l 

that is, we obtain the formula 

p+l 2n 

£>• Sym s (-l) s g p+1 _ s = - E P = 0, 1, • • • , (4.22) 

s=0 a=l 

while this sum vanishes for p — —1. 

From (|4.20p . ()4.21j) . we obtain for R Pji with p > n: 

Rn+l+k,i = -J2q r E ^2n.- S ,2n-l-fc+r(2n - S - i - Sym^^! 
r=0 s=l+i 
2n fc 

+ E ^ +fc ~ l - E qr(-l) t+k - r+n Sym n+fc _ r _,(n - 2 - * + r - 1), 

a=l r=0 

fc = 0,...,n-2. (4.23) 

Two last terms in (J4.23)) pertain to a "regular" part of formulas (J4.12j) and (J4.13|) . i.e., to 
their parts that do not contain T p ^ with p > n — 1. It is straightforward to see that these 
regular parts matches. The only nontrivial thing to do is to check matching "nonregular" 
parts. To make presentation shorter, we use ellipses to denote omitted "regular" parts. 
We have for k = 0, . . . , n — 2: 

n— 2 n—k k—3 

E = E ' ' ' E Rn+l+p,n+l+p-k 

i=0 i=0 p=0 

fc-3 / p n 

= (•••) + E I _ E 9r E T2n-s,2n-l-p+r X 
p=0 \ r=0 s=n+2+p— k 

x (2n -s-(n + l+p-k)- l)(-l)-«+*>+* Sym^^^ 

= (•••)- E T„_ 2+fe _ Sj 2n-l-r E( 2fc -P- S - 4 )(-) S ~ PS y m S -p^-0 
0<r<s<fc-2 \P=r- / 

= (...)- E r n _ 2+fc _ s , 2n _ 1 _ r [[2fc - 2s - A]5 r>s + (- E /C + 2^r, s ) ) 

0<r<s<fc-2 V V a=l / / 

2n 

= (•••)+ E T n - 2 +k-s,2n-l-r E 
0<r<s<fc-2 Q=l 

- E 2[n + k-r-2]T n . 2+k . rt2n . 1 . r , (4.24) 

0<r<fc-2 
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which exactly gives the "nonregular" part in expression (J4.13)) . The lemma is therefore 
proved. 

5 Conclusion 

The obtained answer ()4.9|) has a nice geometrical interpretation in terms of G-functions 
and isomonodromic deformations (see [TOj, JS], CH])- Worth mentioning is the relation 
between this answer and the answer in planar limit, which is a superpotential. These 
geometrical relations deserve special studying. We hope to provide also a geometrical 
explanation for combinatorial identities of type (|4.8jl . 

Worth mentioning is what will be the answer for Fx if we assume the "old" matrix 
model approach with equated chemical potentials on the intervals of eigenvalue distribu- 
tion. It is easy to see that the corresponding answer has form ()4.9|) with the only difference 
that instead of integrating over A-cycles in the matrix Qij, one should integrate over the 
.B-cycles. In the particular case of two-cut solution, this reproduces the Akemann answer. 

Another question pertains to the ^-dependence of F x . The loop equation by definition 
may fix only the part of F± that depends on the potential; in principle, there can be a part 
that depends only on the occupation numbers and which cannot be attained by using the 
Virasoro invariance of the matrix model integral; however, this part must be saturated by 
the planar limit of the matrix model integral. It is nevertheless interesting and important 
to calculate S'-derivatives of the answer 

The author thanks B. Dubrovin, T. Grava, D. Korotkin, A. Marshakov, A. Mironov, 
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